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Abstract 

Integrable Hamiltonians for higher spin periodic XXZ chains are constructed in terms 
of the spin generators; explicit examples for spins up to | are given. Relations between 
Hamiltonians for some t/g(sZ 2 )-symmetric and 17(l)-symmetric nniversal r-matrices are 
studied; their properties are investigated. A certain modification of the higher spin 
periodic chain Hamiltonian is shown to be an integrable t/g(s/ 2 )-symmetric Hamiltonian 
for an open chain. 

PACS: 75.10.Jm, 02.20.Uw, 02.30.Ik 


1 Introduction 


XX Z spin chains have numerous connections with two-dimensional statistical physics and 
(l-l-l)-dimensional quantum field theory. They describe interaction of g-deformed spins 
sitting at the nodes of a one-dimensional lattice. The spin generators S~ and obey 
the commutation relations of the quantum Lie algebra Uq{sl 2 ) [0 


[ 5 +, 5 -] 


sin(275^) 
sin 7 


[S^,S^] =±S^. 


( 1 ) 


We will consider only integrable XXZ spin models. The simplest example in this class 
is a spin-^ chain with the Hamiltonian given by 


H = + (cos7) . 

n 


Higher spin XXZ chains have also been studied @ |> i but only in the spin-1 case an 
explicit expression for the corresponding Hamiltonian has been found Q . 

In the present paper we will construct local integrable Hamiltonians (in terms of the 
spin generators) for higher spin XXZ chains. In Section we recall some facts about 
Gq(s/ 2 )-symmetric universal r-matrix r{X). In Section ^ we construct a 17g(sl2)-sy™™6ti'ic 
local Hamiltonian Hn,n+i- Its properties and properties of the corresponding closed chain 
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Hamiltonian Ti are discussed in Section ^ In particular, we observe that decomposes 

into a ?7(l)-symmetric bulk Hamiltonian Hn^n+i plus a universal local boundary term. Sec¬ 
tion H contains explicit expressions for Hn,n+i for spins 1 and In Section ^ we find 
a family of universal r-matrices and Hamiltonians corresponding to an alternative choice of 
the co-multiplication. In Section ^ we construct another family of (C/(l)-symmetric) uni¬ 
versal r-matrices and local Hamiltonians which contains the Hamiltonian Hn,n+i and the 
corresponding reflection-symmetric universal r-matrix ro(A). In Section ^ we employ our 
construction to establish the integrability of certain (C/q(s/ 2 )-symmetric) Hamiltonians for 
an open chain. Appendix A contains technical details related to r(A) and ro(A). Appendix B 
provides some details on computation of the spin-1 and spin-| Hamiltonians presented in 
Section |^. Appendix C explains a g-trace formula for Hn^n+i used in Section 

For compactness of notations, we will use in the text both the deformation parameter 7 
introduced in and q = . We assume that q is either real or takes values on the unit 

circle. In the latter case q is assumed to be generic, i.e., it is not a root of unity. 


2 Universal r-matrix 

The starting point of the quantum inverse scattering method is the exchange relation 

R{X) L{\ -I- /u) (g) L{pi) = L{n) (g) L(A -|- fi) R{X ), (2) 

with (g) understood as the tensor product with respect to an auxiliary space V (below it is 
the space of 2x2 matrices) and the usual product of operators in the quantum space Sj. The 
i?-matrix belongs to H (g) H. Thus (^) is an equation in H (g) H (g) 

The following L-operator, i.e., an element of V <Si is consistent with the algebra (|I|) 


L(A) = 


sinh[ 7 (A -|- i5^)] i sin 7 e'>''^ S 


sin 7 \ i sin 7 6 S~^ sinh[ 7 (A — iS"^)] 

in the sense that it satishes eq. (0) provided that R{X) = PvR{X), where 


.R(A) = ie'*'^ cr’’" (g) cr -|-ie '^^cr (g) cj’''— 


sm7 


sinh 


17 / 

-I- -^(1 (g) 1 -I- cr 


3 ^^3 


(3) 


(4) 


Here cj“ denote the Pauli matrices, ci^ = cr^ ± iu^, and Py = ^(1 (g) 1 -|- Yla ® 
permutation matrix in V (g) P. 

The L-operator (^) decomposes into two A-independent Borel components. 


(2 sin 7 ) L(A) = — e '^'^L_ , 


(5) 


which can be utilized to define the co-multiplication (a linear homomorphism A : Uq{sl 2 ) —> 
[Uq{sl 2 )]^‘^) in the matrix form, 

A(L±) = L±®L± , ( 6 ) 

where A acts on the quantum space and ® denotes the tensor product with respect to Sj 
and the usual matrix product in V. In explicit form (^) reads 


A(5^) = 0 q-^' + q^' ® , A(S^) = 53 ® 1 + 1 ® 53 


(7) 


Observe that can be obtained by evaluating the quantum space of L(A -|- ^) in the 
fundamental representation, where = q^^a^ and = icr^. In fact, R(X) and L(A) are, 
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respectively, V and V ® representations of a more general object, R(A) e [C/g(sZ 2 )]®^, 
that satisfies the Yang-Baxter equation in (g) ^ (g) 

Ri2(A) Ri 3(A + p) R23(/r) = R23(Ai) Ri3(A + p) Ri2(A) . (8) 

Relation (I) is then a particular case of @ with the first and second spaces evaluated 
in the fundamental representation. Evaluating @ in the V 0 Sj 0 representation (and 
interchanging A with //), we obtain 

r(A) L(A + = L(p)(gL(A + p) r(A), (9) 

with r(A) = R(A)P, where P is the permutation in (g) For a given L-operator, can 
be regarded as a defining equation on the element r(A) E \Uq{sl 2 )]®‘^ which we will call the 
universal r-matrix. Using (^), equation (^) can be equivalently rewritten as follows 

r(A) L±®L± = L±®L± r(A), (10) 

r(A) {e^^L+®L- + e-T'^L_®L+) = {e-'^^L+®L- + e^^L-®L+) r(A). (11) 

In view of (^), the first line implies that for any element ^ E Uq{sl 2 ) we have 

r(A)A(0 = A(0r(A). (12) 


Recall that, for a generic g, the tensor product of two irreducible highest weight Uq{sl 2 )- 
representations of spin S is completely reducible and decomposes into the sum 

25 

Ds®Ds = '^Dj, (13) 

i=o 


where each subspace Dj is a highest weight 17q(sZ2)-niodule with respect to the action of the 
operators A(5=*=) and A(S'^). 

Eq. (p^) implies that r(A) is a function of an operator J such that 


J |j, rn) = j |j, m) (14) 

for any vector |j, m) from Dj. In other words, 

25 

= ( 15 ) 

3=0 


where Vj is the projector onto Dj, i.e., Vk\j,m) = Sjk \j,rn). Taking (p^) into account, one 
can solve equation (|ll]) explicitly ( |^, 0 and see also Appendix A), 


2S j 

r(A) = 7^0 + £ (n 


sin[ 7 (/c — iX)] 

^ sin[7(^ + ^A)] 


Vi 


(16) 


Two obvious consequences of this formula are 


r(A) r(—A) = 1 (g) 1, r(A) r(//) = r(^) r(A). 


( 17 ) 


Introduce a g-analogue of the gamma function satisfying the relation 

{q--q-^)Tq{x) = {q-q-^)Tq{x + l) 
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and normalized such that rq(l) = 1. If |g| ^ 1, this equation can be solved in terms of a 
convergent infinite product. For instance, for \q\ < 1, the solution is given by 


n=0 


1-q 


2n+2 


_ Q2n-\-2x 


(18) 


In terms of the g-gamma function eq. (HD can be rewritten as follows 

^q{J + 1 — iX) rg(l + iX) 


r(A) = 


Xq{J + 1 + *A) rg(i — i\) 


(19) 


with J defined by (p^. Eqs. and (|^ are g-deformations of their XXX counterparts [] 
found in Q] (see also p). 


3 Hamiltonian 


Having a solution R(A) to equation (^) such that R(0) = P (equivalently, r(0) = 1(8) 1), one 
can construct an integrable Hamiltonian for a closed chain in the following way I 0 i 
(the normalization is chosen for later convenience) 


n = i 


. smy 


N 


d 


N 


N 


2 ^Pn,n+l ^^Rn,n+l(A) 
’ n=l 


Hn,n+1 — i 


.sin 7 d 


2j dX 


In r 


n,n+l 


(A) 


A=0 


A=0 


n=l 






n-\-l,n ? 


n=l 


( 20 ) 


with = Rac,! (the periodic boundary conditions). Although Hn^n+i / Rn+i,n) we 

will show below that for the total Hamiltonian we have 


N 


N 


W = ^ Hn+l,n = ^ Rn, 


n+1 


( 21 ) 


ri.=l 


n=l 


The Hamiltonian 7i commutes with the higher quantum integrals of motion which are con¬ 
structed as higher derivatives of Inr(A). Moreover, if the reference state u) for an L-operator 
in question is such that is an eigenvector of r(A), then the corresponding Bethe vectors 

are eigenvectors of the Hamiltonian |M, IT, Q, In particular, this is the case for the 
L-operator ( 0 ), for which the reference state oj is just the highest weight vector. 


Combining (IT) with (^), we obtain a compact formula for the Hamiltonian, 




n,n+l 


Sin 7 

7 


{^q{Jn,n+l + 1) ” ^q{l)) , 


( 22 ) 


where 'l'g(x) stands for the logarithmic derivative of rq(x). 

In order to find the Hamiltonian explicitly in terms of the spin generators we first sub¬ 
stitute eq. (16) into eq. ( 20 ) and derive 




n,n+l 


2S j , 

cosjk 




i=l k=l 


■ U I 
sin 7 /c / 


(23) 


^Actually, [g and ^ dealt with the XXX analogue of R(A). In the XXZ case, P and R(A) do not 
commnte with A and therefore cannot be represented as combinations of the projectors Vj. 


4 



















Next, we have to construct the projectors Vj explicitly. Recall that the Casimir operator of 
the algebra ([^ is given by 

+ + (24) 

sin 7 

Its value in the highest weight representation of spin S is 


sin 75 sin 7 ( 6 '+ 1) cos 7 — cos 7 ( 26 '+ 1) 

^ 2 o ^ 2 * 

Sin 7 2 sin 7 

Applying the co-multiplication (^) to the Casimir operator we obtain an operator that 
acts in the tensor product Ds ® Dg, 

X 5 = ^AC = S+) ® [S- q-^") + S-) ® (5+ (25) 

“*" 4 sin^ ((1 < 8 ) 1 + q^^^ ® cos 7 — (1 (g) + q"^^^ (g) 1 ) cos{'y{2S + 1))^ . 

Here we employed the commutation relations (||) and used (as reflected in the subscript 
of X 5 ) that on Dg ® Dg we have C'(g)l = l(g)C' = (l(g) 1 )C' 5 . It is now easy to see that 


X 5 


sinyj sin 7 (J + 1 ) 
2 sin^ 7 


(26) 


with J defined by (14). Indeed, by construction, X 5 commutes with A{S^). Bearing this in 
mind, d^) follows from computing the action of X 5 on highest weight vectors. 


For a generic q, the eigenvalues of (^) on different subspaces Dj in (^) do not coincide. 
Therefore, we can utilize (26) to construct the projectors Vj by means of the Lagrange 
interpolation. 


25 


p..=n 


2Xg - [l]g[l + l]g 




(27) 




Here we used the g-numbers defined as [k]q = {q^ — q~^)/{q — q~^) = ^ - ^o^e that the 

denominator in (| 2 ^ can be written differently with the help of the identity [j — ^q[j+l + ^]q = 
b 1 qb + + l]i?- We have chosen the first expression to indicate singularities that 

can occur if g is a root of unity. 

Combining (^) with (^), we obtain a local ^integrable Hamiltonian for the XXZ spin 
chain: 


25 r j 

Hn,n+i = (sin 7 ) (^ 7 ] 


j=l k=l 


COS ‘jk 

sin 7/0 


25 

n 

z=o 
17 ^ j 


2 (sin 7 )^ X 5 — sin 'yl sin 'y{l + 1 ) 
sin7(j — 1 ) sin7(j + / + 1) 


(28) 


with X 5 G ijn < 8 > Sjn+i given by (^) which can also be rewritten as follows 

Xs = (1S+S-+, + 1S-S7, + siii7S;sin7S7i ““T;gT+” 

_ +co.75;c°5 7^7i ““AST^‘’ ) ■ (29) 

^We call Hn,n+i local because it is a lattice analogue of the Hamiltonian density. But Hn,n+i is also 
local in the sense that it involves spines only at two nearest sites of the lattice. We hope that this mixed 
terminology will not lead to a confusion. 
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We remark that 2 X 5 is a g-deformation of square of the sum of two spins in the sense that 
in the 7 —> 0 limit eq. (^) simplifies to 

XO = S{S + 1) + S^ SO+i = ^(S° + S0+i)2 , 

where S'® are the generators of 5 / 2 - In this limit eq. (28) turns into the integrable XXX 
Hamiltonian constructed in Q . 


Let us consider the large spin S asymptotics of the Hamiltonian (M)-(23) in terms of 


the spin operator J defined by (W). Recall first that for the logarithmic derivative of the 
(non-deformed) gamma function we have 

If 

+ 1) = ^'(1) + X] + 9) + ■ 


fc=i 


Therefore the large spin approximation of the XXX Hamiltonian is 

Hn,n+l = ln(J + 2) + ■ 


In the XXZ case, if q is real, eq. ( p^ ) (or its counterpart for |g| > 1) yields ^q{x) 
for large x. More precisely, we infer from (p^-(23) that 


(30) 

x\l\ 


^q{j + 1 ) = ^g(l) + ItI ^tanh(/c| 7 |) = \-i\j + + 0{e , 


k=l 


(31) 


where is a 7 -dependent constant. Hence the large spin asymptotics of the XXZ Hamil¬ 
tonian is 

Hn,n+i = J sinh I 7 I -|- const. (32) 

Here we should remind that the spin operator in (^) is not the same as in (|30D, Indeed, it is 
always given by J = jVj but the projectors are different for different 7 . Notice also that 


the correction to the leading order in (M) decays very fast (because tanh(x) = 1 
So, if lylS" is not too small, (|^ is a good approximation even for not too large S. 

Thus, the large spin asymptotics of Hn,n+i in the XXZ case differs from that in the 
XXX case. Moreover, for real 7 , i.e., when |g| = 1, the function 'Lg(x) is not monotonous 
(as seen from (|^)) and does not have an asymptotics at all. 


4 Properties 

Global symmetry: Recall that the co-associativity property, (A (g) 1)A = (1 (g) A)A, leads 
to a natural notion of an n-th power of the co-multiplication: = (A (g) l®’^)Ai”'i 

with A^^i = A. Then the global spin generators for a chain with N nodes are naturally 
introduced ^ as 5^ = (S^) and 5® = Ai^“^i(S'®); which in explicit form reads 


N 


N 




s! 


‘^"+1... g 


(33) 


n=l n=l 

By construction (see (|l2|) and (^) ), the local Hamiltonian Hn,n+i is C/g(sZ 2 )-symmetric, i.e., 
it commutes with A(^) for any ^ G Uq{sl 2 )- Furthermore, it is easy to check that 

[Hn,n+i,S^]=0, for n = 1,... ,iV, (34) 

[Hn,n+i,S^] =0, for n = 1,...,A-1. (35) 
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Since i does not satisfy (|^, the total Hamiltonian enjoys only the C/(l)-symmetry, 

[H,53]=0. (36) 


Actually, the higher quantum integrals of motion also commute with S^. Therefore in the 
presence of a constant magnetic held h the corresponding Hamiltonian, Tih = Ti + hS^, 
remains integrable and C/(l)-symmetric. 

C and P symmetries: Recall that P denotes the permutation in ( 8 ) Since the co¬ 
multiplication A does not commute with P, neither does r’(A), X 5 or Hn^n+i- However, we 
observe that (here and below dependence on 7 is shown explicitly only if it is affected by a 
transformation) 


PX5(7)P = X5(-7), Pr(A,7)P = r(A,-7). (37) 

The hrst relation is obvious from ( p^ and yields the second one upon noticing that both the 
projectors 'Pj(x, 7 ) and the coefficients rj{\,'y) in (^) are even functions in 7 . Further, the 
second relation in (|^) implies readily that 


Pn,n-|-1 Pn,n+l{l) Pn,n-|-1 — P-n,n+l{ T) ■ 


(38) 


Thus Hn,n+i does not have the P (rehection) symmetry. But, if q is real, it has the 
C-symmetry (invariance with respect to the complex conjugation). If jg'l = 1, then eq. (^) 
shows that Hn^n+i has no C- or P-symmetry separately but it has the CP-symmetry. 

Local bulk and boundary terms: As we will see below, Hn^n+i decomposes into two parts 
(which we will refer to as the local bulk term and the local boundary term): 

Hn,n +1 = Hn,n+1 + ” <5^+i) • (39) 

The local bulk term, Hn,n+i, has the following properties (see Section]^ 

Pn,n+l77n,n+lPn,n+l — Hn,n+l j Pn,n+l{p') — Pn,n+l{ T) • (“^9) 


Thus, Hn,n+i is P- and C-symmetric for real q as well as for l^l = 1 . In fact (see Section |^, 
Hn^n+i is a local Hamiltonian associated with the L-operator 


L(A) 


1 

sin 7 


sinh[ 7 (A + iS"^)] 
i S~^ sin 7 


i S sin 7 
sinh[ 7 (A — iS"^)] 


(41) 


Notice that Pn.n+i is not Pg(s/ 2 )-symmetric but only P(l)-symmetric. 

The local boundary term in (^) has the same form for all positive (half-) integer spins. 
In the total Hamiltonian of a closed chain all these boundary terms mutually cancel, hence 


N 


N 


77 = ^ Hn,n+1 = ^2 


n+1 


(42) 


n=l 


n=l 


This together with (|40| ) explains why the two sums in (21) coincide. 

Properties with respect to *-operation: The ^-structure on Uq{sl 2 ) corresponding to the 
compact real form Uq{su{2)) is defined as an anti-automorphism such that 


(5±)* = , {s-^y = s 


(43) 
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with some real rj. Of course, the prefactor can be eliminated by rescaling the generators. 
However it may be convenient to keep it. For instance, we saw in Section ^ that, in the 
spin-i case, it is natural to put 5^ = q^ 2 a^. Then r/ = 1 if |g| = 1 but r] = q if q is real. In 
fact, the choice of rj is not important for our purposes since Hn^n+i and Hn,n+i contain 
only in homogeneous combinations like ( 8 > S~). 

The action of the *-operation extends on a tensor product as ( 8 ) C)* = C* ® C*- H 
must be remarked that properties of the co-multiplication with respect to the action of the 
♦-operation depend on the choice of q, namely we infer from (0) and (H) that f\ 

(A(0)* = A(r) ifczGM but (A(0)* = PA(r)P if k| = 1 • 

Therefore the local Hamiltonian (|2^ has the following properties: 

)* = Hn,n+1 for g G M (44) 

{Hn,n+l) — ^n,n+lHn,n+l^n,n+l fo^ l^l ~ 1 • 

Nevertheless, we see from eq. that the eigenvalues of Hn,n+i are real in the both cases. 
Furthermore, it follows from (|39|)-(p0|) and (D-il) that 

{Hn,n+l)* = Hn,n+l (46) 

for the both regimes of q. In combination with (|^) this implies that the relation 

n* = n (47) 


holds also in the both regimes of q. 

It is worth emphasizing that, for |g| = 1, objects which have the global C/g(sZ 2 )-symmetry, 
like Hn,n+i for n ^ N, are in general not self-conjugate with respect to the ^-operation (H). 
Indeed, if O* = O and it commutes with the global spin generators 5^ given by (p^), then 
it must also commute with (5^)*. This imposes a strong extra condition on the structure 
of O since for |g| = 1 the conjugate of 5^ belong to Uq-i{sl 2 ) rather than to Uq{sl 2 )- Vice 
versa, objects that are self-conjugate for |g| = 1 are in general not C/g(sZ 2 )-symmetric, for 
instance Hn^n+i and TL. 

T symmetry: The definition (^ ® Q* = ^* ® Q* is consistent with the property of the 
matrix transposition, <Si CY = ® C*! if ^ and C are regarded as matrices. Thus the 

♦-operation can be realized as the matrix hermitian conjugation, ^* = (bar denotes the 
complex conjugation). 

From 


and (|38[) we deduce that {Hn,n+i)* = Hn,n+i in the both regimes of q. 
Therefore hnite dimensional matrix representations of the considered above Hamiltonians 
are symmetric matrices, 

{Hn,n+lY = Hn,n+U (#n,n+l)‘ = ^n,n+l , {Hf = H . (48) 


Here the first equality leads to the second and the third due to 


and (|4^, respectively. 


5 Examples (S = 1, |) 

For the spin S' = ^, the Hamiltonian (p^) is very simple: 


Hn,n+i = (cosy) Pi = Xi , 

__ 2 

>i<-structure satisfying (A(5))* = A(^*) for |g| = 1 is the non-compact real form Uq{sl{2,K)), see 


(49) 






( 50 ) 


and the spin generators are given by and = ha^. Using the relations 




a = a e 






= e (7 


Ata'^ 


= COS t + icj^ sin t , 


it is easy to check that (^) acquires the following form 

Hn,n +1 = \{S+S-^, + 5 - 5 ++ 1 ) + (cosy) SlsAi + ^ cosy + 1^(5^ - S^+i). (51) 

The bulk term here is the well-known XXZ deformation of the Heisenberg spin chain. 

For S = 1 the Hamiltonian (^) looks as follows 

Hn,n+1 = (cosy) Pi + ^ „ ((cosy-b 4cos^ y) Xi - (Xi)^) . (52) 

sm2y 4(cosyj'*V / 

In this case the spin generators are given by (|^); they are related to the spin-1 sZ 2 -generators 
as = (cosy) 2 S'^, = Sq. Rewriting Xi and its square as polynomials in the spin 

generators (see Appendix B), we obtain 

Hn,n+l = (f^l - (^^l)' + Fi) + i^ (53 - Sl^^), (53) 

where 

= ^-^(<5;[<5n+i + <5n*5;[+i) + (2cosy- 1 ) 5353 +,, (54) 

Zi cos j 

Fi = 2cosy(cosy- 1 ) 5353 +,-F2(cosy- 1 )^( 5353 +,)2 (55) 

- 2(sin y)2 ((5)))^ (53+,)^) -h 4 2 cos 2y . 


The bulk term of the Hamiltonian (|^ ) (i.e., (|^ ) without the last term) coincides (up to 
a constant) with the Fateev-Zamolodchikov (FZ) Hamiltonian Q. The FZ Hamiltonian is 
usually written in a slightly different way in terms of the spin-1 sZ 2 -generators; we have 
chosen the above form since it allows for better comparison with the 5 = | case. 

In the limit y —> 0, (K)-(^) simplifies to the well-know spin-1 XXX Hamiltonian: 


trO _ 5 vO 1 /'y0'i 2 1 cO C' 

^n,n+l — 4 4 1^17 — 4 ^ 


n+1 - i (s: 


. cO 12 I 3 

n 2 ' 


In the 5 = 1 case the Hamiltonian (^) is given by 


ffn,n+i = (cosy) Pi 


sin 3y 
sin 2y 


P2 + 



(56) 


Rewriting the corresponding expression (^) in the polynomial (with respect to the spin 
generators) form (see Appendix B), we obtain 


H, 


1 


n,n+l — 


12(cos y)3(l -|- 2 cos 2y)3 
1 


^12(cosy) (ri3)3 -|- (5cos4y — cos2y -|- 2) (fis) 




I -|- 2 cos 2y 


2 cos y 2 


(fis Q + Q 113) + ) + i^(53 - 53 +,), (57) 


64 cos^ y 2 


.smy 


where Q and F 3 are real symmetric polynomials in 53 and 53+, (see Appendix B) and 


^n+i + S-S+^,) + (1 + 2cos2y)|^ 5353+,. 

2 .Z O COS I 


(58) 
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Let us remark that in both Hamiltonians (53) and (57) the two “most non-diagonal” 
terms are powers of an operator Li which is quadratic in generators. To clarify this fact, we 
observe that, like for 5 = ^, the corresponding Casimir operators (U) are actually quadratic: 


-(5+5- + 5 - 5 +) + (cos 7 )(5^)2 = 2 cos 7 for 5 = 1, 

i(5+5- + 5 - 5 +) + (cos7)2(5=^)2 = ^( 17 cos 2 7-2) for 5 = ^ 


(59) 

(60) 


(the last expression differs from Cs by a constant). Eq. (|5^ is due to the first relation 
in (1^ ; eq. (^) follows from the relation (|107D and the identity 2(sin 75 ^) ^ = 1 — cos 275 ^. 


6 Another co-multiplication 

As we saw in Section]^, the co-multiplication operation A plays a key role in the construction 
of the Hamiltonian Ti. However, eq. (|^ does not exhaust possible definitions of A for Uq{sl 2 )- 
In fact, a homomorphism Aq, : Uq{sl 2 ) \Uq{sl 2 )]®‘^ such that 

A„(5+) = 5++ g(+“+^)^'®5+, A„(53) = 5^^1 + 1053 (61) 

satisfies all the properties of the co-multiplication for any real a. We had before a = 0 which 
is the most “symmetric” choice. But we may prefer a = 1 (or a = —1) if we want the regime 
q < 0 to be on equal footing with g > 0 in the sense that the terms q(+"+r)s3 (|^) take 

only real values. 

Notice that Aq is obtained from the co-multiplication (@) by twisting: 

A„(0 = A(e) (E„)-1 , (62) 


This twist is rather specific in that it preserves the co-associativity of the co-multiplication 
(more general twists give rise to the so-called quasi-Hopf algebras |H|; see ||T^ for their 
applications in integrable spin models). 

In order to construct an integrable Hamiltonian using the new co-multiplication A„ in 
the same way as we used A in Sections ^ and we should first find an L-operator, I/q(A), 
that satisfies (1^ with Aq. For this purpose we observe that 


Fa = (</>a ® </>a) ^ A((/)a) where (l)a = G ^ ■ (63) 

Therefore we can rewrite the l.h.s. of (P) for La{X) as 

Aa{La{X)) = {(j)a ® La{X) (pa^) i4>a <8) </>a) 

which makes it obvious that La{X) satisfies (^) with A„ if it is related to the L-operator (^ 
as follows 


LaiX) 


(j)jL{X)(j)a 

1 / sinh[ 7 (A -|- i5^)] 

sin 7 I i sin 7 5+ 


i sin 7 e'^^g 2 g“'S'^ 5 \ 

sinh[ 7 (A — i5^)] j 


(64) 

(65) 


Since the map ^ 4'a^^4‘a is an automorphism of Uq{sl 2 ), it is clear that the exchange 

relation (^) holds for La{X) with the same i7-matrix as for L(A). Consequently, the Bethe 
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ansatz equations for La{X) coincide with those for L{X) (the reference state w has also not 
changed). 

Substituting (|^ into (^), we find a universal r-matrix for the new L-operator: 

fo(A) = {(l)a < 8 ) r{X) {c/)a < 8 ) </>«) = Fa r{X) {Fa)~^. ( 66 ) 

The second equality is due to relation (^) and the property (p^). The same twists relate 
the corresponding solutions of the Yang-Baxter equation (^), i.e., 

Ro(A) = ((/)„ (g) 4>a)~^'R-{X){(t>a ® 4>a) = Fa R(A)(Fq,)“^ . (67) 


Notice that the first equality is consistent with (^) since (j)a in the fundamental representa¬ 
tion is just a constant. Further, evaluating the r.h.s. of ( |^ in Y (g) representation, we see 
that La{X) can be constructed also as a twist by 2 x 2 matrix, 

La{X) = faL{X)f-\ = 

Let us underline that existence of the two ways of constructing La{X) and, as a consequence, 
of the relation 

[L{X),cl)afa\=Q 

is due to the property of the universal r-matrix (^) applied to ^ = <pa (indeed, q'^4’afa is 
A{(j)a) evaluated in V (gjij). 

According to (|^), the local Hamiltonian corresponding to the universal r-matrix ( | 66 [ ) is 
^S+l = « C+l)"' ^n,n+l (C <+l) = ■ (68) 


This transformation does not modify the S' = ^ Hamiltonian (49) since in this case cj)a is 
trivial. But already for S = 1 we find (see Appendix B) 


u(«) _ TT _ 

“n,n+l J^n,n+l — 


1 — COS (ay) 
2 cos 7 


{T„,„+i,S;)S;)+a-i 


. sin(a 7 ) 


2 cos 7 


(69) 


where [,] and {,} stand for commutator and anticommutator, respectively, and T^^n+i = 
-|- S“S^_^;^). Notice that the terms on the r.h.s. of (| 6 ^) are non-diagonal. 

Finally, the total Hamiltonian is given by 


N 


w'“’ = E'^2+i = (®= 


\-l 


TC^a = Fa FL {Fa 


\-l 




n=l 


N 

H' 

n=l 


Fa^ q 


The <h-twist here follows easily from (^). The .T-twist yields the same result because FL 
commutes with (36) and hence with ^aFa = ■ Since 77^“^ and FC are related by a 

twist they have the same set of eigenvalues; this agrees with the fact that the Bethe ansatz 
equations have not changed. 

Let us conclude this section with a remark: eqs. (|^), (|^ and (|^ may appear to suggest 
that Lq(A) and ra{X) are related to another twisted co-multiplication, A 0 (^) = 6~^ A{^)6, 
where 9 = (j)a0 (pa- But fails to satisfy the necessary property of a co-multiplication []T|, 
(e (g) l)A(^) = where e is the co-unit. So Ag is not a co-multiplication of a (quasi) Hopf 
algebra. (Moreover, A^ is not co-associative.) 
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7 From Hn,n+i to 

Consider now another family of L-operators, 

•p _ 1 / sinh[ 7 (A + i5^)] i siny 5“ 

^ sin 7 y i sin 7 5+ sinh[ 7 (A — i5^)] 

which are obtained from the L-operator (|3|) as twists by certain 2 x 2 matrices, 

L0{X)=K^^L{X)Kx, = 


(70) 


(71) 


For P = j this gives the L-operator (^) (which is most often used in applications of the 
Bethe ansatz). The map is an automorphism of Uq{sl 2 ) but, unlike 

the case treated in Section P, it is A-dependent. Therefore, the L-matrix corresponding to 
Lfj{X) differs from (Q). Namely, as seen from ([7T|), exchange relation fl) holds for Lp(X) with 
Rp{X) = (1 ( 8 > K^^)R{X){Kx < 8 > 1), i.e., RpiX) = PRj 3 {X), where 


Rp{X) = i fj+ ® a- + i 


a 


) (j^ + • 


sin 7 


• sinh 


+^(i' 


i 1 -|- cr^ (g) 


(72) 


In order to find a universal r-matrix for L,g(A), we will apply the approach which we used 
in Section Namely, we observe that tp\ = G is a complementary twist to (0) in 

the sense that 

[L(A),iFAV^A] = 0 

as follows from the property of the universal r-matrix ( 0 ) for C = (fx (in fact, Kx is just (px 
evaluated in the fundamental representation; hence Kxp\ is A( 73 a) evaluated in 

Thus, instead of the twist (0D in the auxiliary space V, the L-operator (|70|) can be 
obtained as a twist in the quantum space Sj, 


Lp{X) = pxL{X) 


G ^ . 


(73) 


Substituting (^) into (^) and taking again into account that [r(A), A((^^)] = 0, we find 
a universal r-matrix for L^(A), 

f^(A) = (1® <^a)?’(A) (<^a^ 


1 ). 


(74) 


Unlike eq. (pbD this relation is not a twist. Notice however that the corresponding solutions 
of the Yang-Baxter equation (0) are related by a twist, 


R-/3(A) = (1 ® <y9A)R(A)(l ® 




(75) 


It must be stressed now that L^(A) does not possess a decomposition of the type dD-d), 
and rp{X) does not commute with A (or Aq,), i.e., (lY) does not hold for r, 3 (A) for a generic 
As a consequence, f(g(A) does not have a representation of the type (|l5|). Moreover, for /3 7 ^ 0 
we have in general [r( 3 (A), 773 (^)] 7 ^ 0 (except for the fundamental representation in the case 
/3 = 7 ). Nevertheless, the general recipe for constructing a local integrable Hamiltonian 
applies (because the Yang-Baxter equation for R,g(A) is valid). So we substitute (74) into 
the formula (^) and derive 


_ ZJ _ 4 

^n,n+l ~ ^n,n+l 1 


: P sin 7 ^^3 


27 




(76) 
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Thus the new local Hamiltonian differs from only in the local boundary term. Hence 

the total Hamiltonian H = Yln=i n+i coincides with H. Note that the Bethe ansatz 


equations, describing the spectrum of 7i, have also not changed, which is not entirely trivial 
since the corresponding i?-matrix has changed. The reason is that in the derivation of 
the Bethe ansatz equations the non-diagonal entries of (72) appear only in the so-called 
“unwanted terms” that cancel each other 11 , 0 , 1 - 

In Section ^ we asserted that the local Hamiltonian Hn,n+i decomposes into a local 
boundary term and a local bulk term Hn,n+i which is associated with the L-operator (^). 
Now this is obvious from ( |70|) and (|7^ if we put /3 = 7 . Since we have found the corre¬ 
sponding universal r-matrix, we are in a position to prove the properties of Hn^n+i stated in 
Section |^. For brevity we denote ro(A) = f 7 (A). 

As seen from (^) for /3 = 7 , the auxiliary ii-matrix associated with the L-operator (^) 
is P-symmetric, i.e., it commutes with the 4x4 permutation matrix Pv- The corresponding 
universal r-matrix, ro(A), has analogous properties. Namely, (as we prove in Appendix A) 
ro(A) satisfies the following relations 


Pro(A)P = ro(A), ro(A, 7 ) = ro(A,- 7 ), (77) 

(?’o(A))‘= ro(A), ro(A)ro(-A) = 1(8) 1, (78) 


where P is the permutation in ( 8 ) and t denotes transposition. The last equality follows 
from the formula (0), the first relation in (H) for r(A), and the relation [r(A), A((y 9 ;^)] = 0. 
Taking logarithmic derivative of (0) at A = 0, we establish the symmetries (|40| ) of the local 
bulk term Hn,n+i- The first relation in (|7^ ) yields in the same way the second relation 
iniD. 


8 Open chain 


As we saw in Section 77^,1 is the only term in Ti which does not commute with the global 
spin generators S^. Omitting it, we obtain a Hamiltonian for an open spin chain, 


Af-l 


Af-1 


^n,n+l = - Sn) 


(79) 


n=l 


n=l 


which is apparently tAj(s/ 2 )-symmetric, i.e., [7i',S^] = [7{',S^] = 0. It is however not 
immediately evident whether the Hamiltonian (|7^ remains integrable. 

Let us refer to the sum on the r.h.s. of (|7^ as the bulk Hamiltonian, Ti'. The remaining 
part can be called the surface term; it is a sum of {N — 1) local boundary terms we dealt 
with before. As we discussed above, the bulk Hamiltonian Ti' corresponds to the L-operator 


(41) and it is only t/(l)-symmetric. The fact that adding the surface term to 7i' restores 
the t/g(s/ 2 )-symmetry was observed for spin ^ in ||l^ and for spin 1 in 0 ]- Integrability 
of the corresponding total Hamiltonians was established in [18, in the framework of 

boundary integrable lattice models. We will prove below that 7i' is integrable for the higher 
spins as well. 

Let us briefly recall the construction of an integrable Hamiltonian for a chain with bound¬ 
aries IQ. Let R(A) be a solution of the Yang-Baxter equation (^), and T{X) be a monodromy 
matrix obeying the exchange relation (§) with P-matrix PR,(A). Introduce a boundary mon¬ 
odromy matrix, Z{X) = T{X)}C~ {X){T{—X))~^. The boundary matrix,^ (A) e Sj, repre- 

“^Strictly speaking, IC~ (X) is an element of (g) but with trivial second component. We had a similar 
situation in Section g where the twist tfi\ had trivial first component, see eq. ([t^. 
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sents non-periodic boundary conditions (R and JC~ are close analogues of the bulk and 
boundary scattering matrices |^). Now, if /C“(A) satisfies the so-called reflection equation, 
namely 

R(A - //) (/C-(A) ® 1) (R(-A - ^))-i (1 0 JC-{fj,)) (80) 

= (1 ® /C-(/i)) R(A + m) (/C-(A) ® 1) (R(-A + /r))-^, 


then Z{X) also satisfies this equation. Using this fact one can show that a special trace of 
Z(X), r(A) = [X)Z[X), is a generating function for quantum integrals of motion if the 

boundary matrix /C'''(A) G Sj satisfies a “dual” reflection equation |T^, ^|, 


R(-A + ((/C+)‘(A) ® 1) R(-A -/I -26) (10 (/C+)*(^)) (81) 

= (10 (/C+)^(fi)) R(-X -fi-26) ((/C+)*(A) ® 1) R(-A + fi) , 

and R has the following properties 

PR(A)P = R(A), R(A)R(-A) = /9i(A)(1® 1), (R(A))* = R(A), (82) 

(R(A))'i (R(-A - 26))*^ = p2(X) (10 1), (83) 


where 5 is a constant, and /9i(A) and / 92 (A) are scalar function. 

With all these conditions, an integrable Hamiltonian is given |T^ by the following ana¬ 
logue of the formula (|^) (we keep the same normalization as in (p0|)) 


iV-l 




n=l 


47 dX 


A=0 


+ 


trQ(/C(]~(0) hp^jv) 
tr/C+(0) 


• sin 7 d ~ 

Pn,nH-l 

Z7 CLA 


A=0 


(84) 


Deriving (|^ one assumes that R„^n+i(0) = P and /C“(0) = 1, which is consistent with (^). 

Let us try to identify the Hamiltonian (|79|) as a particular case of (pi]). First, we can 
put hn^n+i = Rn,n+i if we choose R(A) = R.y(A), where R.y(A) = ro(A)P is given by (^) 
with /3 = 7 . For this R-matrix the properties ( p^ follow from ([7^; in particular, the second 
relation (unitarity) holds with /9i(A) = 1. The crossing unitarity (| 8 ^) holds for R.y(A) with 
<5 = i. 


The derivation of the reflection equation (|80D for )C~ does not use the conditions 
(82). So, let us look first for an R-matrix for which the reflection equation has a trivial 
solution, IC~(X) = 1. In this case (p^ turns into 


R(A - p) (R(-A - p))-^ = R(A + /u) (R(-A + /r))"^. (85) 


A solution to this equation is given by R(A) = r(A)P, where r(A) is the universal r-matrix 
we discussed in Section |^. Indeed, ( p5| ) follows from the second relation in (0)- Using this 
observation, we can find a solution to the reflection equation for any i?-matrix R/ 3 (A) given 
by ([ 7 ^). Indeed, substituting (|75D in (p5[), we derive 

R;3(A - p)( 10 ifl) (R/3(-A - p))~^ (1 0 ipl^) (86) 

= (10 V’-,,) + /i) (1 0 V>1) (R/3(-A + . 


Multiplying this relation by A(p^‘^) = ((p\0(p\) ^ (which commutes with R/ 3 (A)), we bring 
it to the form of (30). Thus, a solution of the reflection equation for R /3 is 
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Notice that for /3 = 7 the i2-matrix R is symmetric with respect to 7 (the second relation 
in ([77|)). Therefore in this case we have another solution (which can be derived directly from 
(| 8 ^) by applying the permutation and then multiplying by ® 

/C-(A) = . (87) 


This is the boundary matrix we need since its derivative in (U) gives exactly the S\ term 
in d^. 

In general, solutions of the reflection equation and of its “dual” are independent of each 
other. However, as was noted already in [^], there exist several isomorphisms that allow us 
to construct /C+ if we know JC~. In particular, it is easy to see that a possible solution for 
(11) is /C+(A) = (/C-(-A - 6))K For /C' 


given by (^) this yields 
/C+(A) = . 


( 88 ) 


It turns out that substitution of /C^(0) = into (|^) gives exactly the term in (^). 

To prove this assertion, we first observe (see Appendix C) that for the q-ii&ce of the local 
Hamiltonian Hn^n+i we have 



77)1,n+l) — PS ln+1 ) 


(89) 


where ps is a scalar constant. 

Taking into account the relation (^^ between Hn,n+i and Hn^n+i-, we infer from (|^) 
that trQ(( 7 ^'^o Hq^at) = (tr ( 7 ^'^^) i^^^ 5^ + const. Replacing now q with q~^ and using that 
Hn^n+i is an even function of 7 (|40|), we find that 


tro(/C+(0)Ho,Ar) _ trg(g ^^0 _ .sin7 3 

tr/C+(0) “ trg-2S3 “ ' 2 ^ 


(90) 


holds (up to an additive constant). Thus, the boundary matrix (^) gives the Sfj term 
in d?^). This completes the proof that the [/q(s/ 2 )-symmetric open chain Hamiltonian (^) 
is integrable for all (half-) integer spins. 


Conclusion 


In summary, we have constructed explicitly (in terms of the spin generators) higher spin 
closed chain Hamiltonians for two families of XXZ-type L-operators including the two 
L-operators most often used in the literature. We have investigated properties of these Hamil¬ 
tonians, described their interrelations, and discussed the connection with 17g(s/2)-symmetric 
open chain Hamiltonians. 

We have emphasized a key role of the underlying quantum algebraic structure, especially 
the universal r-matrix and the co-multiplication, for constructing integrable Hamiltonians 
and investigating their symmetries. The technique presented in this paper can be applied 
also for constructing the higher quantum integrals of motion. 

Our construction does not in general extend to the case when q is a root of unity. However 
it appears that we could allow = 1 if p is sufficiently large in comparison with the spin S. 


For instance, examining the denominator of (^) we see that the projectors Vj do not become 
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singular at roots of unity if p > 85. If the ratio p/S is not large, then an appropriate 
modification of the presented construction is needed; it will be certainly of interest and use. 
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Appendix A 


Here we give some technical details on the universal r-matrices used in the text. First we 
recall the derivation of the universal r-matrix (jl^) along the lines of (^, ^ (see also 1 . 11 ). 
The explicit form of the off-diagonal entries of equation (0) is 

r(A) (e^T'^5± ® ® 5±) = ® 0 5±) r(A). (91) 

A solution to this equation is unique up to a scalar factor |^. Notice that the co-multiplication 
dD satishes the following relations 

[5" ® q^', A(5")] = [g-^' ® 5", A(5")] = 0 , e = ± . 

Therefore, for the highest/lowest weight vectors we deduce that 

\j,±j) = g±{j)\j + l,±{j + l)), g"^"(8)5=^ |j,±j) = h±{j)\j + l,±{j + l)) , (92) 


where g±{j) and h±{j) are scalar functions. Furthermore, we have 


(5± ® g^" + g±(2+2Hg-S3 ^ 5 ±) 1 ^-^ ^ ^ q-S^ ( 93 ) 

+^±(2+2i)(^53^^±^2S3^^_253^^-2S3^ |j, ±j) = (1 ® g^^^) A(5±) |j, ±j) = 0 . 

Applying now (|^ to \j,j) (for the upper signs) or to \ j, — j) (for the lower signs) and using 
(jl^), (p^) and ( 1 ^, we deduce that 

rj+i(A) (e^T-^ - g±(2+2i)eT7A) = _ q±{2+2j). (94) 


Both relations in ( [^ ) yield the same functional equation on rj{X) (which arises also for 
universal r-matrices in the lattice sine-Gordon model [^| and in the lattice Virasoro alge¬ 
bra [p5|) 


rj+iW = 


sin[ 7 (j -M - iX)] 


KA) 


(95) 


sin[ 7 (j + 1 + iA)] ^ 

Upon imposing the normalization condition r(0) = 1, we obtain the expression (p^. A proof 
that this r(A) does satisfy all the relations in (|Tl] ) is given in i 0 - Of course, one is still 
free to multiply r(A) by a scalar function, p(A), such that p{0) = 1. For instance, the spin-^ 
representation of R{X) given by (§) corresponds to p{X) = sin[ 7 (l -|- iX)]/ sin 7 . 

Consider now the universal r-matrix ro(A) introduced in Section It is a solution to 
equation (|^ for the L-operator (^). Unlike the previous case, this L-operator does not have 
a Borel decomposition of the type (^). Therefore the off-diagonal entries of (^) give us four 
A-dependent relations: 


ro(A) (5± ® q^^ + ® 5±) = (e±^^5± ® g^' + g”^' ® 5±) ro(A), (96) 

ro(A) (5^ (g) q~^^+ q^^ ® 5^) = (g) 5^ -|- 5^ (g) g“‘^^) ro(A). (97) 
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It can be proven along the lines of that solution to this equation is unique up to a 
scalar factor. As we have already shown in Section |^, this solution is given by ro(A) = 
(1 (g) (g) 1), where r(A) solves (^. 

Now we observe that ri(A) = Pro(A)P and r 2 (A) = (ro(A))* solve the same set of equa¬ 
tions By the above mentioned uniqueness, this implies that ri{X) = Cj(A)ro(A), 

where Cj(A), i = 1,2 are scalar functions. Since P^ = 1 (g) 1, and {{roYY = vq, we conclude 
that Ci(A) = ±1. Imposing the condition ro(0) = 1 (g) 1, we have to put Ci(A) = 1. Thus, we 
have proven the relations on the l.h.s of (|7^ ) and (|7^). Employing the first of them and the 
property (^) of r(A), we prove the second relation in (^) as follows: 


ro(A, -7) = (1 ® e-^^^')r(A, 0 1) = (1 ® e-T'^^')Pr(A,7)P(eT'^^' ® 1) 

= ® l)r(A,7)(l ® eT'^^")P = P(1 ® ® 1)P 

= Pro(A,7)P = ro(A,7). 

In the second line we used that [r(A), A(e'^'^‘^^)] = 0. 


Appendix B 

Here we provide some details on computation of the Hamiltonians in the cases of spin 1 and 
spin |. For 5 = 1, a matrix representation of the spin generators is (as was discussed in 
Section |^, one has a freedom of rescaling 5^ —> with any real 7 7 ^ 0) 

/0a0\ /100\ 

S+ = 0 0 a , S- = {S+Y^ 53 = 0 0 0 , (98) 

\000/ \00-l/ 


where t denotes transposition and a = ^J2 cos 7 . Since {S^Y — function of is a 

polynomial in of a degree not exceeding two. In particular, we have 


sm(tS'^) = S'^ sint, cos{tS'^) = 1 — 2(5’'^)"^ sim - 


(99) 


Denote T = i(S'+S'^_|_j^ + 5'^ 5'^+i)- With the help of formulae (p^ we rewrite d^) as follows 

Xi = -t-(cosy) i-I-2 (cos 7)^^ , (100) 

A = ^-2(cos7)' + (cos7)'(5353+i) + (sm7)'((S^53+,)2-2(53)2-2(53+1)2) 

+ i(sin 27)(53 - 53+1 + ^(53)253+^ - ^53(53+1)2) . 


Substituting (|100D into (p^ , applying several times formulae 
Hamiltonian in the following form 


and (M), we obtain the 


H, 


1 


n,n+l 


1 '-'^2 1 


(T)2-(T A + A T) 


4 cos 7 V cos^ 7 ' ' cos 7 
-1-4 + 2 cos 27 -|- (cos 27) (5353+1 — (53s'3+i)2) 

- 2(sin7)2((53)2 + (S3+i)2)) + i!^ (53 - 5, 


( 101 ) 


^+ 1 ) 
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Deriving the first term here we used the following analogue of relation (50) 


( 102 ) 


Next, we observe that the following identity holds (as can be checked directly in terms of 
matrices) 


with Fn as in 


and 


O(T„)-i = 0-T + T0+, 

(103) 

2+“ 0037 ) 5353+1 . 

(104) 


Substituting (103)-(104) with a = 0 into (101), we obtain the Hamiltonian (p^). 

The Hamiltonian discussed in Section ^ is obtained from (|101j) by the twist (|^). 

Notice that the term (T)^ in (1^) is not affected due to (1^). The diagonal part of 
Hamiltonian apparently commutes with the twist. So the only part of (1^) which changes 
is the second term. It transforms according to (1^), which yields the Hamiltonian 


Consider now the case S' = |. The spin generators are given by 


S+ = 


/ 0 ai 0 0 \ 

0 0 02 0 

0 0 0 oi 

V 0 0 0 0 / 


/ 


s- = {s+Y 


s-^ = 


0 0 


0 


1 


0 0 -i 


0 \ 
0 
0 


(105) 


V 0 0 0 -I / 


with oi = (2 cos 27 + 1)2 and 02 = 2 cos 7 . Any function of is a polynomial in S^ of 
degree not exceeding three. In particular. 


sin(tS3) = (2sin^ + isin3l)S3-l(sin3i)(S3)3, 
cos(tS^) = cos ^ + - sin ^ sint — (sin - sint) (S^)^ . 

The Hamiltonian (p^) for S = | takes the following form in terms of X 3 


H, 


72 , 71+1 — 


■^4(cos7) (Xs)^ 


(106) 

(107) 

(108) 


4(cos 7)^(1 + 2 cos 27)3 

— (Xs)^ (13 + 20 cos 27 + 8 cos 47 + 2 cos 67 ) 

2 

+ X 3 (68 cos 7 + 48 cos 87 + 23 cos 67 + 7 cos 7'y + cos 
2 

This Hamiltonian can be rewritten in the same way as we did for S = 1 using, in particular 
relations 


(|106D and (|107|) , and appropriate analogues of ( |102D and (|103|) . The final form 
is given by (|5^, with 


Q = 


1 

1 ^ 


(1 ( 8 ) 1) (815 — 5832 cos 7 + 2352 cos 27 — 3888 cos 87 


+2542 cos 47 — 1620 cos 67 + 1600 cos 67 — 324 cos 7'y + 467 cos 87 ) 
sin 'S'n'S'n+i (106 + 37cos 7 + 186 cos 27 + 29 cos 87 
+88 cos 47 + 13 cos 57 + 28 cos 67 + 5 cos 77 ) 
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+(sin (28 — 53cos 7 + 54cos 27 — 31 cos 37 

+28 cos 47 — 11 cos 57 + 10 cos 67 — cos 7j) 

+ ^(sm|)^ ((5'^)^ + (*S’^+i)^) (412 + 1687cos7+ 606cos27 

+ 1085 cos 37 + 268 cos 47 + 433 cos 67 + 58 cos 67 + 83 cos 7'y), 

2 

COS 'y 

Fs = —-—(1 (g) 1) (10328 + 17080 cos 27 + 9071 cos 47 + 3220 cos 67 + 621 cos 87) 

-^(sin 5'^5’^+i (8364 + 5020 cos 7 + 12752 cos 2j + 2499 cos 37 + 7150 cos 47 
9 2 

+ 1023 cos 57 + 3320 cos 67 + 447 cos 77 + 814 cos 87 + 83 cos 97) 

--(sin (5'^'S'^+i)^ (536 — 472 cos 7 + 864cos 27 — 347cos 37 + 650cos 47 

— 163 cos 57 + 272 cos 67 — 65 cos 77 + 78 cos 87 — 9 cos 97) 

— — (sin ('S'n'S'n+i)^ (50 — 200 cos 7 + 48 cos 27 — 168 cos 87 + 45 cos 47 
9 2 

— 76 cos 57 + 28 cos 67 — 20 cos 77 + 13 cos 87) 

-(sin27)2 (^(5^)2 + (53+1)2) (169 + 274cos27 + 106cos47 + 27cos67) 

+-(sin ^)2 ((53)353+1 + 53(53+1)2) (444 — 2516cos 7 + 560cos 27 — 2037cos 87 
+334 cos 47 — 1137 cos 57 + 200 cos 67 — 417 cos 77 + 46 cos 87 — 85 cos 97). 


The equivalence of (105) and (57) as 16x16 matrices has been verified with the help of the 
program Mathematica. 


Appendix C 


Here we explain the origin of the equation (^) that was important for our discussion on the 
open chain Hamiltonian in Section ^ 

Recall that, by definition (see e.g., ||T^, ^]), a Hopf algebra A possesses the antipode map 
s : A ^ A which is an anti-homomorphism consistent with the co-multiplication and the 
co-unit in the sense that m((s (g) fd)A(^)) = m((i(i (g) s)A(^)) = e(^) • 1 (here m : A®"^ —> ^ is 
the multiplication). Assume that there exists an element x G .4. which realizes the square of 
the antipode (which is a homomorphism) as an inner automorphism, i.e., s(s(^)) = X^X~^ 
for any G .4. Then, as was proven in pffl. 


(109) 


tri((x (g)l)6) 

belongs to the center of A if an element b G A®'^ satisfies [b, A(^)] = 0 for any G 
For A = Uq{sl 2 ) the antipode consistent with the co-multiplication (^) is given by 

s(5±) = 

It is easy to see that in this case x = ■ Since the universal r-matrix r(A) has the property 

(H), we can apply (|109D and infer that its g-trace, trj^((( 72 'S'^ 0 l) r(A)), belongs to the center 
of Uq{sl 2 )- Consequently, the same holds for the local Hamiltonian Hn,n+i constructed as 
in (p0|). Being evaluated in an irreducible representation, the g-trace of Hn^n+i becomes just 


a constant, as was stated in eq. (^) 
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